function Fig3_a()
Fcolor='k';
Bcolor='w';
cMap=getColor('Magma');

lable2 = {'-5/6','-19/24','-3/4','-17/24','-2/3','-3/5','-7/12','-13/24','-1/2','-11/24','-5/12','-2/5','-1/3','-7/24','-1/4','-5/24','-1/6'};
nu2 = -[5/6, 19/24, 3/4, 17/24, 2/3, 3/5, 7/12, 13/24, 1/2, 11/24, 5/12, 2/5, 1/3, 7/24, 1/4, 5/24, 1/6];
ind2 = [61,79,62,78,110,112,64,76,114,75,66,113,111,72,68,71,69];

result = num2cell(zeros(1,length(ind2)));
parfor i = 1:length(ind2)
    result{i} = CS1_0_extraction_2D_Phi(ind2(i),0);
end

figure();
for i=1:length(ind2)
    [nr_sub,nc_sub] = size(result{i}.phi);
    tmpy = meshgrid(1:nr_sub)';
    tmpx = meshgrid(1:nc_sub);
    tmpy = reshape(tmpy,[],1);
    tmpx = reshape(tmpx,[],1);
    [fitobj,~,~]=fit([tmpy,tmpx],reshape(result{i}.phi,[],1),'poly11');
    map_bkg = meshgrid(1:nr_sub)'*fitobj.p10 + meshgrid(1:nc_sub)*fitobj.p01 ;
    map_bkg = map_bkg - mean(mean(map_bkg));


    phi = (result{i}.phi - result{i}.avg - map_bkg)*1e3;
    win = gausswin(100,100*2.5/4)*gausswin(100,100*2.5/4)'; %std=4nm
    phi_plot = conv2(phi,win,'same')/sum(sum(win));
    
    subplot(1,length(ind2),i)
    imagesc(result{i}.rx(1,:),result{i}.ry(:,1),phi_plot)
    ax = gca;
    title(lable2(i))
    ax.YDir='normal';ax.XDir='normal';
    axis equal;axis square;
    ax.XTick=[];
    ax.YTick=[];
    ax.Color=Bcolor;
    ax.XColor=Fcolor;
    ax.YColor=Fcolor;
    colormap(ax,cMap);
    ax.CLim=[-0.5,0.4];

end

end